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Abstract. We study the integrability of polynomial vector fields using Galois theory of 
linear differential equations when the associated foliations is reduced to a Riccati type folia- 
tion. In particular we obtain integrability results for some families of quadratic vector fields, 
Lienard equations and equations related with special functions such as Hypergeomctric and 
Hcun ones. The Poincare problem for some families is also approached. 



Introduction 

Given a polynomial differential system in C^, 

dx . ^. . dy . , . 

— = x = P{x,y), — = y = Q{x,y}, (1) 

with P,Q E C[x, y], we consider its associated differential vector field 

x = Pix,y)^ + Qix,y)§^, (2) 

whose integral curves are intimately related to the solutions of system ([1]). These solutions, 
taken as curves on the plane and leaving for a while its time- dependence, constitute its so- 
called foliation and satisfy the first order differential equation 

^ dx P{x,y)' 

This expression is often written as a Pfaff equation 

= 0, (4) 

where Q = Q{x, y)dx — P{x, y)dy is the corresponding differential 1-form. The connection 
between integral curves of the vector field X and solutions of ^2 = is clear: 

• geometrically, it is given by Lx^ = 0, which means that the vector field X is tangent 
to the leaves of the foliation (the orbits) defined by (j4j); 
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• dynamically, the general solution of equation (jl]), H{x,y) = ctt, is given by a first 
integral H of the original vector field X, that is, a non-constant scalar function which 
remains constant along any of its solutions {x(t),y(t)). Since Lx^ = 0, this is equiv- 
alent to say that X{H) = and to the existence of a suitable scalar function / such 
that n = fdH. 

From a geometrical point of view, we focus our attention on invariant algebraic curves, that is, 
polynomial integral curves of the vector field X. This is the natural framework where Darboux 
Theory can be applied. Dynamically speaking, we restrict ourselves to first integrals H{x, y) 
which are Liouvillian, i.e., written as a combination of algebraic functions, quadratures and 
exponential of quadratures in C{x,y), the field of rational functions on x,y. As we will see 
later, Galois Theory provides very useful and powerful tools to approach it. 
Two classical problems remain still open for complex polynomial fields: 

(i) Concerning the existence of invariant algebraic curves of system ([T]) or, equivalently, 

of algebraic solutions of the foliation equation 
(ii) About the existence of Liouvillian first integrals for systems ([T]) (or, in other words, 
to determine when the general solution of equation (|3]) is Liouvillian). 

For general polynomial vector fields, problems (i) and (ii) are very difficult and we are still 
far from obtaining an effective method to decide whether a given arbitrary polynomial field 
has or not an invariant curve or admits a Liouvillian first integral. In fact, problem (i) is 
connected with the (also unsolved) classical Poincare problem, which seeks for a bound of 
the degree of the invariant algebraic curves as a function of the degree of the vector field (or 
of the associated foliation defined by (|3])). It is known that Darboux Theory and adjacent 
results, as the ones due to Prelle-Singer and Singer ED] , provide connections between this 
two problems. 

In this work we are concerned with the study of the (Darboux, Galois) integrability of 
some families of equations of type (jl]) inside the complex analytic category, that is, when the 
original vector field X defining f2 is complex polynomial or can be reduced to it. Precisely, 
we restrict ourselves to those systems which can be reduced to a Riccati type equation 



ai,a2 and 03 being rational functions with complex coefficients. For Riccati equations there 
is a very nice theory of integrability in the context of the Galois theory of its associated 
second order linear differential equation. This is, in our opinion, a natural framework where 
several known results concerning integrability of Riccati equations (|5]) should be considered. 
Kovacic provided in 1986 (see |30]) an effective algorithm which allows to decide whether an 
equation(|5]) has got an algebraic solution or not. And, additionally, a theorem of Liouville 
(see [21]) proved that the existence of an algebraic solution is exactly the definition of the 
integrability for (JS]) in the context of the Galois theory for linear differential equations. Thus, 
for foliations of type ^ problems [i] and [ii) are equivalent and Kovacic algorithm becomes 
an extremely powerful tool to approach them. 

In some sense, this work can be considered as a very particular case of the Malgrange 
approach to the Galois theory of codimension-1 foliations [391 SOI E], that is, for Riccati 
codimension-1 foliations on the complex plane. Our target is not to obtain general theoretical 
classification results, but several effective criteria of integrability for such foliations. We 
will provide integrability criteria for some families of polynomial quadratic vector fields and 
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some Lienard equations involving special functions, allowing this to recover previous results 
established by several authors. For instance, we will solve completely the integrability problem 
for the family of Lienard type equations of the form (see [IB]) 

yy' = (a(2m + k)x^'' + b{2m - k)x'^~^^^)y - {a^mx^^ + cx^^ + h'^m)x^'^"^^~^ , 

with a, 6, c, m, k complex parameters. It is known that when this equation comes from a 
polynomial vector field the constants m and k must be rational numbers. 

We also approach the Poincare problem for some particular families of systems (see Propo- 
sition [XI] and Theorem 13. 4p . 

The paper is structured as follows. To make it as self-contained as possible, we introduce 
in Section [1] a basic background about Galois theory of linear differential equations and 
Darboux theory of integrability of polynomial vector fields. In Section [2] we remind some 
useful properties concerning Riccati equations and Section [3] is devoted to applications. For 
completeness we include two Appendixes about Kovacic algorithm and some special functions. 

1. Two NOTIONS OF INTEGRABILITY FOR PLANAR POLYNOMIAL VECTOR FIELDS 

1.1. Darboux theory of Integrability. We give a very brief overview of Darboux's inte- 
grability ideas [20], his terminology and some essential results. 

Let us consider a vector field ([2]) and an irreducible polynomial / G C[x,?/]. The curve 
/ = is called an invariant algebraic curve of vector field ([2]) if it satisfies 

/l/=o = 0. 

This condition is equivalent to the existence of a polynomial K G C[x,y], called cofactor, 
such that 

X(/(x, y)) = P{x, y)^ + Q{x, y)^ = K{x, y)f{x, y), 
or, equivalently, that 

^ = X(log(/)) = K. (6) 

From this expression it follows that the curve / = is formed by leaves and critical points of 
the vector field X = [P, Q) defined by ([2]). If the polynomial system ([1]) has degree d, that is 
d = max{degP, degQ}, then we have that degK < d — 1, independently of the degree of the 
curve f{x,y) = 0. From definition (E]) it follows that if the cofactor K vanishes identically 
then the polynomial / is a first integral of the vector field X. In terms of the associated 
foliation, this invariant curve / = is a particular solution oi y' = Q/ P and Qdx — Pdy = 0. 

An analytic C-valued non-constant function fi is called an integrating factor of system ([1]) 
if the expression 

^ = X(log(/i)) = -V-X, 

holds, where V ■ X = (dP/dx) + (dQ/dy) is the divergence of the vector field X = {P,Q). 
In case that the domain of definition of X is simply connected, from the integrating factor fi 
it follows that 

Hix, y) = J y)P{x, y) dy + ip{x), 
is a first integral of X, provided that OH /dy = —fiQ. 
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To ensure the existence of a first integral for a system ([T]) is, in general, a very difficult 
problem. In [20], Darboux introduced a method to detect and construct first integrals using 
invariant algebraic curves. Namely, he proved that any planar polynomial differential system 
of degree d having, at least, d{d + l)/2 invariant algebraic curves, admits a first integral or 
an integrating factor which can be obtained from them (see also Jouanolou [27] for a study in 
a general context for codimension-1 foliations). Furthermore, Darboux's original ideas have 
been improved by taking into account the multiplicity of the invariant algebraic curves (see 
[T9] for more details). Related to them some other invariant objects have been introduced 
(see [H]). They are the so-called exponential factors: given h,g E C[x,y] relatively prime, 
the function F = exp (g/h) is called an exponential factor of the polynomial system ([T]) if 
there exists a polynomial K G C[x,?/] (also called cofactor) that satisfies the equation 

^(^) -X(^). K. (7) 



F \h, 

It is known that if h is not a constant polynomial then /i = is an invariant algebraic curve 
of ([1]) of cofactor Kh satisfying that X{g) = gKh + hK. 

The following theorem (coming originally from Darboux) shows how the construction of 
first integrals and integrating factors of ([2]) can be carried out from its invariant algebraic 
curves. 

Theorem 1.1. Let consider a planar polynomial system ([T]) of degree m, having 

• p invariant algebraic curves fi = with cofactors Ki, for i = 1, . . . ,p and 

• q exponential factors Fj = exp{gj/hj) with cofactors Kj, j = 1, . . . ,q. 
Then the following assertions hold: 

(a) There exist constants Xi, Xj G C not all vanishing such that 

i=i j=i 
if and only if the multivalued function 

f^^...f^-F}^...F}^, (8) 

is a (Darboux) first integral of system ([1]). 
(6) There exist constants Xi, Xj G C not all vanishing such that 

p q 



XiKi + ^ XjKj + divX = 0, 



i=i j=i 

if and only if the function defined by is a (Darboux) integrating factor of X . 

For more recent versions of Theorem 11.11 see [371 ESj and for some generalizations see [3^ [7] . 

Functions of the form ([8]) are called Darboux functions. We say that the polynomial system 
([1]) is Darboux integrable if it admits a first integral or an integrating factor which is given by 
a Darboux function. 

Remark 1.2. Prelle and Singer [18] showed that if system ([1]) admits an elementary first 
integral then it admits an integrating factor which is the n-th root of a rational function (a 
slightly improved version of this result can be found in [33| Corollary 4]j. Later, Singer in [50] 
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showed that if system ([T]) admits a Liouvillian first integral then it has an integrating factor 
which is given by a Darhoux function. This is an important argument to motivate sentences 
like "Darhoux functions capture Liouvillian integrahility" or 'Liouvillian first integrals are 
either Darhoux first integrals or integrals coming from a Darhoux integrating factor" . 

Given a polynomial system ([1]) of degree m, the computation of all its invariant algebraic 
curves becomes a complicated problem since nothing is known a priori about the maximum 
degree of these curves. This makes necessary to impose additional conditions either on the 
structure of the system ([T]) or on the nature of such curves (see for instance, [ISl [HI [101 [H] 
or references therein). This difficulty has motivated the study of different types of inverse 
prohlems of the Darboux theory of integrability [121 [HJ [13 [IS, [HI [3S]- 

We finish this subsection with a remark about a geometrical meaning of the integrating 
factor, pretty known to people coming from the Sophus Lie mathematical community. It is 
an established fact in Fluid Dynamics that integrating factors arise as a density in stationary 
planar regimes: the equation V ■ {^iX) = is the continuity equation for the field of velocities 
X, with density function yU = /i(a;, y) (considered in the context of the Symplectic geometry 
and Hamiltonian dynamics). 

Remark 1.3. Let fi = fi{x,y) he an integrating factor of the vector field defined in some 
domain of the plane. Then the vector field X is a hamiltonian vector field with respect to the 
symplectic form Q = fidx A dy (this form degenerates only at the zeros of fi). In fact from 
ix^ = dH , we ohtain 



which is hamiltonian with Hamilton function the first integral H . It is straightforward to 
verify that the symplectic form Vt is invariant under the action of the flow of X , i.e., 



From this point of view, the dynamics of the vector fields in the plane can he formally con- 
sidered as an Ergodic theory prohlem: the existence of an invariant measure, the one defined 
hy the associated integrating factor of the flow. 

1.2. Picard-Vessiot theory. Picard-Vessiot theory is the Galois theory of linear differential 
equations. We will just remind here some of its main definitions and results but we refer the 
reader to [19] for a wide theoretical background. 

We start recalling some basic notions on algebraic groups and, afterwards, Picard-Vessiot 
theory will be introduced. 

An algebraic group of matrices 2 x 2 is a subgroup G C GL(2, C) defined by means of 
algebraic equations in its matrix elements and in the inverse of its determinant. That is, 
there exists a set of polynomials Pj G C[a;i, . . . , 2:5], for z G /, such that A G GL(2, C) given 




Hence, the vector field can he rewritten 




Lx^ = dix^ + ixd^ = dix^ = ddH = 0. 



by 
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belongs to G if and only if Pi (xn, X12, X21, X22, (det = for alH G / and where det A = 

XuX22 — a;2ia;i2- It is said that G is an algebraic manifold endowed with a group structure. 
Recall that a group G is called solvable if and only if there exists a chain of normal subgroups 

e = Go<Gi< . . .<Gn = G, 

satisfying that the quotient Gi/Gj is abelian for all n > z > j > 0. 

It is well known that any algebraic group G has a unique connected normal algebraic 
subgroup G° of finite index. In particular, the identity connected component G^ of G is 
defined as the largest connected algebraic subgroup of G containing the identity. In case that 
G = G^ we say that G is a connected group. Moreover, if G° is solvable we say that G is 
virtually solvable. 

The following result provides the relation between virtual solvability of an algebraic group 
and its structure. 

Theorem 1.4 (Lie-Kolchin) . Let G C GL(2, C) be a virtually solvable group. Then, G^ is 
triangularizable, that is, it is conjugate to a subgroup of upper triangular matrices. 

Now, we briefly introduce Picard-Vessiot Theory. 

First, we say that (/C, ' ) - or, simply, /C - is a differential field if /C is a commutative field 
of characteristic zero, depending on x and ' is a derivation on /C (that is, satisfying that 
(a + b)' = a' + b' and {a ■ b)' = a' ■ b + a ■ b' for all a,b G /C). We denote by C the field of 
constants of K,, defined as C = {c G /C | c' = 0}. 

We will deal with second order linear homogeneous differential equations, that is, equations 
of the form 

y" + b,y' + boy = 0, 61,60 e/C, (9) 
and we will be concerned with the algebraic structure of their solutions. Moreover, along this 
work, we will refer to the current differential field as the smallest one containing the field of 
coefficients of this differential equation. 

Let us suppose that 1/1,1/2 is a basis of solutions of equation ([9]), i.e., yi,y2 are linearly 
independent over K, and every solution is a linear combination over C of these two. Let 
C = /C(|/i,|/2) = /C(|/i, 1/2, i/'i, 1/2) be the differential extension of /C such that C is the field of 
constants for /C and C In this terms, we say that C, the smallest differential field containing 
/C and {1/1,1/2}, is the Picard-Vessiot extension of /C for the differential equation 

The group of all the differential automorphisms of C over /C that commute with the deriva- 
tion ' is called the Galois group of C over /C and is denoted by Gal(£//C). This means, in 
particular, that for any a G Gal(£//C), cr(a') = (cr(a))' for all a G £ and that a{a) = a for 
all a E KL. Thus, if {1/1,1/2} is a fundamental system of solutions of ([9]) and a G Gal(£/A^) 
then {o"?/i, (TI/2} is also a fundamental system. This implies the existence of a non-singular 
constant matrix 

(c rf) eGL(2,C), 

such that 

a [yi y2) = {(r{yi) (y{y2)) = [yi 1/2) A„. 
This fact can be extended in a natural way to a system 

^(yi vA - f ^(^1) ^(2/2)\ _ fyi y2\ 4 
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which leads to a faithful representation Gal(£//C) — t- GL(2,C) and makes possible to con- 
sider Gal(£//C) as a subgroup of GL(2, C) depending (up to conjugacy) on the choice of the 
fundamental system {yi,y2}- 

One of the fundamental results of the Picard-Vessiot Theory is the following theorem 
(see [281 EI]). 

Theorem 1.5. The Galois group Gal(£//C) is an algebraic subgroup o/GL(2, C). 

We say that equation is integrable if the Picard-Vessiot extension £ D /C is obtained 
as a tower of differential fields IC = Cq C Ci C ■ ■ ■ C Cm = C such that Li = £j_i(?7) for 
i = 1, . . . ,m, where either 

(i) r] is algebraic over that is r] satisfies a polynomial equation with coefficients in 

{ii) rj is primitive over that is rj' G A-i- 

(Hi) rj is exponential over that is rj'/r] G 

Usually in terms of Differential Algebra's terminology we say that equation ([9]) is integrable 
if the corresponding Picard-Vessiot extension is Liouvillian. Moreover, the following theorem 
holds. 

Theorem 1.6 (Kolchin). Equation is integrable if and only if Gal(£//C) is virtually 
solvable, that is, its identity component (Gal(£//C))° is solvable. 

For instance, for the case a = in equation (Q, i.e. y" + by = 0, it is very well known [281 
[3T| B9] that Gal(£//C) is an algebraic subgroup of SL(2,C) (remind that A G SL(2, C) 4^ 
A G GL(2, C) and det A = 1). For a more detailed study see the Appendix |X1 

2. Some remarks about Riccati equation 

Ricatti equation is probably one of the most studied equations in Dynamical Systems. 
Its role in the study of the Darboux and Picard-Vessiot integrability leads us to devote 
this section to some of its properties. Even though these results are known, their proofs 
have been included for completeness. We divide these properties in two types: the first one 
(see Subsection 12. ip concerning transformations leading a general second order differential 
equation into a Riccati equation (written in the so-called reduced form). Remind that this 
becomes the starting point of the celebrated Kovacic algorithm (see Appendix . A second 
type, more Darboux-like, that studies first integrals and integrating factors for a Riccati 
equation (Subsection 12. 2p . 

2.1. Transformations related to Riccati equations. It is known that any second order 
differential equation can be led into a general Riccati equation through a classical logarith- 
mic change of variable (see, for instance, [IHl l2Sj)- The following proposition recall it and 
summarises some other related transformations. 

Proposition 2.1. Let IC be a differential field and let consider functions ao(x) , ai(x) , a2{x), 
r{x), p{x), bo{x), bi{x) belonging to K, that, for simplicity, will be denoted without their 
explicit dependence on x. Consider now the following forms associated to any second order 
differential equation (ode) and Riccati equation: 
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(i) Second order ode (in general form) : 

y" + biy' + boy = 0. 
(ii) Second order ode (in reduced form) : 

e = Pi- 

{Hi) Riccati equation (in general form): 

v' = ao + aiv + a2V^, 02 7^ 0. 
{iv) Riccati equation (in reduced form): 

w' = r — w'^, 

Then, there exist transformations T, B, S and TZ leading some of these equations into 
other ones, as showed in the following diagram: 

f ' = ao + aiV + a^v^ w' = r — nP' 



n 



i/" + 6ii/' + 6oy = o— ^r = pe 

The new independent variables are defined by means of 

f a'2 ai\ I I y' 

\2a^ 2a2 / ^2 «2 1/ 

S: y = ^e-^^^^'^^ 7^: ^ = |. 

and the functions r, p, bo and bi are given by 

r = — (ao + aia + a2a^ — a') , 



/ ^2 ai \ 1 
2a2 J 0.2 



2 ' 2- - ' ' ^ ' 



&i = - ( fli + — ) , ^0 = aoa2, 

Proof. The proof is quite standard. 

[T] : Applying the change w = a + fiw we get the equation 

01 + (i'w + (iw = ao + a\oi + a\(iw + a2a^ + 2a2a/3w + a2/3^w^ 

that, regrouping terms, leads to 

w' = — (ao + a\a + a2a^ — a') + + 2a2« — w + a2/3u'^. 
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Since 02 7^ we can take /3 = — 1/02 and, therefore, 02/? = —1. Having this into 
account, the value of a satisfying that the coefficient in w vanishes is given by 

1 fP' 
a = TT ~ '^1 



2a2 \(3 

The expressions for a, f3 and r follow straightforwardly, 

p ' V2«2 2a2 / 02 

Moreover, it is clear that a, /3 and r belong to /C. 



\E\. Imposing a = and taking /3 = — 1/02 in transformation T we have w = —wja^ 
and we obtain the Riccati equation 

/ ( 2 
w = — ao'22 + c^i H ] w — w . 

Performing now the change of variables w = (logy)' (or, equivalently, v = —a2y'/y) 
we obtain the differential equation y" + biy' + b^y = with 

bi = - ( ai + — ] , bo = aoa2. 

Obviously, bo and 61 belong to /C. 
[5]: The change of variable y = fi^, with fi = fi{x) and ^ = 'C(^); l^ad us to 



In order to obtain the equation ^" = we need to impose 

2^ + 60 = 0, ^ + 60- + bi = -p, 
/i 1^1^ 



which gives rise to 



/i = e-^/^°, p = ^ + ^-bi. 



4 2 

Moreover, it is straightforward to check that p G /C. 

[TZ]: This is a particular case of transformation [B] with the particular choice 
ai = and 02 = —1. 

Finally, composing the transformations provided by [B], [R] and [S]: 

-a2V = —, y = ^e"^/^o^ i)^ = -fa^j^^\ ^ = w, 

y V «2/ ^ 

we recover the result given by [7~], 
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which imphes that, in some sense and taking p = r, the diagram commutes: S o B = TZ oT . 

□ 



From this Lemma, it follows that the function v is algebraic over /C if and only if the 
function w is also algebraic over K,. Furthermore, in such case, the degree over K of both 
functions v and w is the same. 

It is known that a Riccati equation f|T2|) has an algebraic solution over /C if and only if the 
differential equation ffTOl) is integrable in a Picard-Vessiot sense. In this situation we say that 
the Riccati equation is integrable over /C. We notice that Kovacic algorithm (see Appendix |A|) 
starts from an equation in form f|TT]) . 

2.2. Integrating factor and first integrals for Riccati vector fields. We briefly show 
some relations between the existence of invariant curves of a certain type of vector fields and 
the integrability, via Kovacic algorithm (see Appendix |A]), of its associated Riccati foliation. 

From Singer |50] (see Remark II. 2p we know that if a planar polynomial vector field ([2]) 
admits a Liouvillian first integral then it has also an integrating factor given by a Darboux 
function. However, very few results are known about the relation between the existence of 
an algebraic invariant curve of a general planar vector field and the Liouvillian integrability 
of its foliation. 

Let us consider a family of planar vector fields of the form 

X = {p{x)-q{x)w^)— + q{x) — , (18) 

with p{x),q{x) G C[x] complex polynomials. Introducing an independent variable t, usually 
called time, we can associate to them the following system of differential equations 

lii = p{x) — q{x)w'^, 
X = q{x), 

where we denote by ' = d/dt. Its foliation, governed by the equation 

, dw p(x) — q(x)w^ p(x) 9 

w = — = -—^ — = — — — — w , 

dx q{x) q{x) 

is a Riccati equation given in reduced form w' = r{x) —w"^ with r = p/q E C(x). Next lemma 
shows that the integrability of this "Riccati foliation" is closely related to the existence of an 
algebraic invariant curve of its vector field (fT8|) . A similar approach for this problem can be 
found in [SUES]. 

Lemma 2.2. Let wi = Wi{x) be a solution of a Riccati equation in reduced form 

w' = r{x) — w"^, 

with r(x) =p{x)/q{x) G C[x]. Then the associated vector field ffTKjl has an integrating factor 
given by 

-2 J wi{x)dx 

Pi{w,x) = j ■ -— ^. (19) 
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Proof. It is straightforward to check that if Wi{x) is a solution of = p/q — w"^ then 
it holds X{fi) = Kifi with fi{w,x) = —w + Wi{x) = and Ki = —q{w + Wi{x)). In 
addition, Fi{x) = e^/'^^W'^^ satisfies X{Fi) = LiFi with Li = -quoi. Notice that X has 
divergence divX = —2qw, and additionally it holds —2Ki + 2Li + divX = 0. Proceeding as 
in Theorem II. ![ vector field (fT8|) admits the integrating factor 

_p2 g-2 / wi{x)dx 

as it was claimed. □ 



Remark 2.3. We stress that the result in Lemma \2.2\ is independent of the nature of the 
solution w = Wi{x). In general, the integral J ui{x)dx is an abelian integral. 

The important fact is that, conversely, Picard-Vessiot theory and in particular, Kovacic 
algorithm, supply information about first integrals and integrating factors of the equation 
w' = r{x) — w"^ from the knowledge of some of its solutions, ^1,^2,^3. Indeed, from the first 
three cases in Kovacic algorithm |30] (the integrable ones) one obtains the following types of 
first integrals (see Weil [51] and Zol§dek [57]). 

Proposition 2.4. The following statements hold. 

Case 1: One has two possibilities: 

- If only Wi G C(x) then X admits a first integral of DarbouxSchwarz-Christoffel 
type. 

- Ifbothwi,W2 € C(x) then X admits a first integral of Darboux type. In particular, 
from Lemma lKE we can construct two integrating factors /ii and /i2 and so ^1/^2 
is a first integral of X . Thus, we have 

Hiw X) = + ^^2(3:)) ^ f\(w2(x)-wr(x))dx] 

^ ' ' {-w + wi{x)) 

Case 2: If Wi is a solution of a quadratic polynomial then X admits a first integral of hyper- 
elliptic type. 

Case 3: If all wi,W2,W3 are algebraic over C{x) then X admits a rational first integral. 

The following result characterises the rational integrability of the polynomial vector fields 
that we are considering in this work. 

Corollary 2.5. The Galois group of f ITT]) is finite if and only if its corresponding planar 
polynomial vector field has a rational first integral. 

Proof. Notice that Gal(£//C) is finite if and only if we fall in case 3 of Kovacic Algorithm or 
in case 1 of the type 

Gal(£//C) = I [| c" = l|. 

Hence, only remains to study this last case (the cyclic one). Let .^1, .^2 be solutions of ^" = p^. 
Then there exists g G C(x) such that ^1 = g^ and ^2 = fl'""- We define Ui = 002 = ^2/^2, 
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and we obtain 

1 ^' 1 ^' 

Ui = , UJ2 = . 

n g n g 



Thus, the corresponding vector field f[TT|) admits the first integral 

H,{u,x) = Z^^i±i^e/(---)'^^ = y, g^. 

-u + ui 1 g 

-cu -\ 

n g 



1 / -ngu - g 



and so also admits the first integral 

g^ \-nguj + g'J 

which completes the proof. □ 



Remark 2.6. Let P G C[x] he a polynomial of odd degree. It is known that the planar 
polynomial vector field x = 1, y = P{x) +y'^ (with associated foliation y' = P{x) +y'^) is not 
integrahle, that is, it has no invariant curves, since it falls in case 4 of Kovacic Algorithm |3Uj . 

2.3. Riccati foliations. Let us recall some well-known geometrical properties of Riccati 
foliations defined by planar polynomials vector fields (see, for instance, [32]). Although we 
are not using these properties along the paper, we include them for completeness. 
Let 

Vt = q{x)dy - {pi{x) + p2{x)y + p:>,{x)y'^)dx, pi,gGC[x] (20) 
be the 1-form defining a Riccati foliation on the complex plane. Since the equation = is 
the projective version of the corresponding second order linear differential equation defined 
over the vector bundle x (i.e., the fibre is projectivized to P^, see subsection 12. ip . 
the Riccati equation f2 = is defined in a natural way over (x,?/) G P^ x P^. 

The singular points of = are the zeros of the polynomial q{x) and, possibly, the 
point at infinity x = oo G P^. Moreover, these singular points are exactly the poles of 
the coefficients of the associated second order linear differential equation. We define d : = 
max{deg{pi),deg{p2),deg{p3),deg{q) — 2). Thus the point x = oo is a singular point if and 
only if deg{q) is less than d + 2. In Kovacic algorithm, which applies to the reduced form 
of the second order linear differential equation (see Appendix this set of singular points 
is denoted by F. Therefore it seems natural to call it in the same way also here, that is, 
r = {xi,...,Xr}. Thus the Riccati foliation is holomorphic on (P^ — F) x P^, because the 
singular points of the associated linear differential equation are exactly the set F, given by the 
singularities of their coefficients. We notice that the ("singular") sets {xj} x P^ are invariant 
by the foliation and are usually called invariant fibres because they are already fibres of the 
fibration tt : P^ x P^ — P^, (x, y) t— )■ x. This fibration is transversal to the Riccati foliation 
since fibres T(x) := 7r^^(x) = {x} x P'^, with x non singular, are global and transversal to 
the foliation, i.e., transversal to all the leaves. Over any of these transversals the holonomy 
group of the foliation is defined representation 

7ri(pi-F,xo)^Diff(T(xo)), 
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(where Diff(T(xo)) is the group of diffeomorphisms on the transversal), given by hfting the 
loops in the fundamental group to the leaves of the foliation, this is, by solving the Riccati 
equation with initial conditions and final points on the transversal T{xq). As the Riccati equa- 
tion is the projectivization of a second order linear differential equation, the holonomy group 
must be the projectivization of the monodromy group of the linear second order equation 
acting on the vector space fibre of the meromorphic vector bundle x C^, {x} x ~ C^. 
By fixing a base of fundamental solutions, this can be considered as the space of solutions of 
the linear differential equation. Hence, as the monodromy group is represented by the linear 
group GL{2, C), the holonomy group is represented by the projective linear group PGL{2, C), 
the Mobius transformations 



The Riccati foliations are the most well-known class of a family of foliations, the projective 
foliations, with holonomy group represented in the projective group. For Riccati foliations, 
the holonomy group is contained in the Galois group of the foliation either in the Malgrange 
approach |39l HHl |9] or in the Lie-Vessiot-Kolchin approach O |5]. In fact, if the singular 
points of the associated linear differential equations are singular regular ones, then the Zariski 
adherence of the holonomy group is the Galois group of the foliation, because in this case the 
Zariski adherence of the monodromy group of the linear differential equation is the Galois 
group of the associated linear differential equation. 

The critical points of the associated vector field, i.e., zeroes of q{x) and of pi{x) +P2{x)y + 
Ps{x)y'^ are obviously contained in the invariant fibres. For general Riccati foliations there 
are two critical points on the invariant fibre for each point in F. For some special Riccati 
foliations there are no critical points. This is the case, for example, of Riccati foliations 
given in reduced form by 1] = q{x)dy — (p(x) -|- q{x)y'^)dx (with p and q relatively primes), 
corresponding to the field (|T8|) . 



In this section we analyse some examples involving integrability and non-integrability of 
some families of Riccati planar vector fields or planar vector fields whose foliation can be 
reduced into a Riccati form. 

3.1. Quadratic polynomials fields. The study of the integrability of the quadratic poly- 
nomial vector field 



with Ojj, hi J G C is, in its general form, a hard problem. One of its possible approaches is the 
so-called linear-quadratic case, when one of the two components is a polynomial of degree 
one. In [36| Prop. 3] it is proved that its study around a finite equilibrium point (the origin) 
can be reduced to consider two families of systems. Using the notation introduced therein, 
we refer to these families as (Sl)-type, 



7ri(P^ - r,a;o) ^ PGL(2,C). 



3 



Applications 



X = a2ox'^ + CLuxy + ao2?/^ + clio^ + a^iy + aoo, 
y = b2oX^ + buxy + bo2y'^ + biox + boiy + boo, 



X = X, 

y = £x + \y + b2ox'^ + biixy + feo2y^ 



(SI) 
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and (S2)-type, 

y = ex + Xy + h2QX^ + huxy + 602^^- 

In [36], the authors prove that the hnear-quadratic systems having a global analytic first 
integral are those satisfying: 

(01) 602 = A = 0. 

(61) 602 = and A = -p/q E Q", 

in the case of (Sl)-type systems and 

(02) ho = fco2 = A = and ebn ^ 0. 
(^2) ^20 = &11 = A = and £602 7^ 0. 
(C2) fell = A = and 620 7^ 0, 

for (S2)-type systems. Furthermore, they also provide the explicit form of the corresponding 
first integrals. It is important to notice that all of them are of Darboux type and, therefore, 
Liouvillian. 

Our aim in this example is to show that these results can be recovered using arguments 
coming from the Galois theory of linear differential equations. We start first with the (Sl)- 
case, whose associated foliation is given by the Riccati equation: 

dy , . /A + 6iixA 602 2 /oiN 

— = (£ + 6202;)+ h/ + — y- (21) 

ax \ X I X 



Kx) = 7-r + ^^ -=^^^i^(^02. + ^(l-A)), /. = ^, (22) 



By Lemma [2.11 this equation can be transformed into the reduced form w' = r{x) — w"^, with 
1 K ^ V - 1 

4 X .yj^— '±^20^02 

provided &11— 4^20^02 7^ 0, and into the form ^" = r{x)C,- This equation is a Whittaker equation 
(see Appendix B) to which one can apply the Martinet-Ramis Theorem (see theorem IB.2p . 
This Theorem asserts that such Whittaker equation is integrable if and only if at least one 
of the following conditions is verified: 

±K ± /i G ^ + N, 

or, equivalently (and more suitable for the expressions derived of n and fi) 

2 (k ± /i) G 2Z + 1. 

In our case one has that 

A/^n - 4620&02 

so for (Sl)-type systems conditions (ai) and (61) read, respectively, 2(k ± /i) = 1 G 2Z + 1 
and 2{k + /i) = (1 + (p/q)) + {—p/q) = 1 G 2Z + 1. Therefore Galois Theory recovers the 
integrability result asserted in [36| Thm.l]. 

Let us consider now a (S2)-type system, namely, 

X = y, 

y = £x + Xy + b2ox'^ + buxy + 6o22/^ 
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with foliation given by the differential equation 

^ = (A + hix) + {ex + 620a;') - + bo2y. (23) 
ax y 

This equation falls in one of the following situations: 
(i) A = 611 = yields to a Bernoulli equation 

= Ux + b2oX } - + 6022/, 
ax ^ ' y 

which corresponds to cases (62) and (02). 
{ii) If e = 620 = we obtain the linear equation (and, of course, integrable in a Liouville 
sense) 

^ = (A + hiix) + ho2y- 
ax 

This possibility is not taken into account by Llibre and Vails [3S] since this equation 

is not, strictly speaking, in Riccati form. 
{iii) &20 = 602 = A = and ehu 7^ (case (02)) gives rise to dy/dx = bux + exy~^ , which is 

a separable equation (and a Bernoulli as well) and whose solutions are all Liouvillian. 
(iv) If 602 = we obtain a Lienard equation, 

^dx ^ ^^^^^ y+ {ex + b2ox'^) , 
that will be considered more deeply in a forthcoming section. 

3.2. Families of orthogonal polynomials. Recall that the Hypergeometric equation, in- 
cluding confluences, is a particular case of the differential equation 

y" + + = 0' ^ e C, L = ao + aix, Q = bo + bix + b2X^ . (24) 

It is well known (see, for example, [13]) that classical orthogonal and Bessel polynomials are 
solutions of equation flMl) for suitable values of a^ , bj and A. Namely, 

• Hermite Hn-, 

• Chebyshev of first kind T„, 

• Chebyshev of second kind f/^, 

• Legendre P„, 

• Laguerre L„, 

• associated Laguerre 

• Gegenbauer 

• Jacobi V^'"^ 

• Bessel 

where 
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Family 


Q 


L 


A 


Hfi 


1 


-2x 


2n 


Tn 


1 - 


—X 




Un 


l-x' 


— 3x 


n{n + 2) 


Pn 


l-x' 


-2x 


n{n + 1) 


Ln 


X 


1 — X 


n 


T (m) 


X 


m + 1 — X 


n 




l-x^ 


— {2m + l)x 


n{n + 2m) 


1 n 


l-x" 


V — m — {m + V + 2)x 


n{n + 1 + m + u) 


Bn 


•> 

X 


2{x + l) 


—n{n + 1) 



Integrability conditions and solutions of equation flMj) can be obtained applying Kovacic 
algorithm (Case 1 of the algorithm). Besides, they can also be achieved via Kimura and 
Martinet-Ramis Theorems and the parabolic cylinder equation (see [221 [23111] )• 

As a consequence, from this table, we obtain the following result. 

Proposition 3.1. We consider Q, L and A as in the previous table. Then, for any /i 7^ 0, 

the planar quadratic polynomial vector field 



dv 
dx 



A 

Q, 



Q + {Q'- L)v + 



(25) 



VJx) 



Q'-L u 2 

V H V . 

Q Q 



where Vn{x) is any orthogonal 



has invariant algebraic curves of the form ftv 
polynomial associated to X,Q,L and n eN. 

Proof. The associated (Riccati) foliation of system (l25l) is written into the form 

dv A 
dx ji 

Performing the change v = fiv equation fj26|) becomes 



(26) 



dv 
dx 



A + 



Q'-L^ 1 



Q Q 

and can be led into the form (124|) through the transformation v = —Qy'/ y (see transformation 
B of Proposition 12. ip . For fixed A, Q and L, let Vn be any orthogonal polynomial solution of 
equation (124|) . Then Vn = — (QV^) / (fiVn) is a rational solution of equation (126|) . Therefore 
the curve f — f „ = is an invariant curve of the vector field (l25|l for any n G N. □ 



According to Proposition 12.11 equation f l26|l can be reduced to the form ^" = with 

(g) "q + (2q) ' 

f — 

and ^ = VnG'' is a solution for any n G N. Additionally, we notice that we fall in Case 1 of 
Kovacic algorithm. 
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3.3. Lienard equation. Let us consider first order differential equations whose associated 
foliation can be expressed into the Lienard form 



yy' = f{x)y + g{x), 



(27) 



with y = y{x) and rational functions f{x) and g{x). We are concerned with the problem of 
obtaining criteria on f{x) and g{x) such that equation fl27j) can be led into a Riccati equation. 

This is a difficult problem and, as far as the authors know, only partial answers have 
been given to it. In what follows we give some examples of such results coming from the 
handbook [IB] and papers [TT| [T^. 

A first example is given by the 5-parametric family (1.3.3.11 in [15] ) 



yy 



(a(2m + k)x^^ + h{2m - k)x 



y 



(a^mx^^ 



cx 



2k 



h m)x 



2m-2fc-l 



(2^ 



a, b, c, m, k being complex parameters. In order that (l28i) come from a polynomial vector 
field we have that m and k must be rational numbers (see [16|). The change w = x'^, 



X 



[z + ax^ + bx ^) leads (EHD into the Riccati form 



{—mz^ + 2abm — c)w'{z) = bk + kzw + akw'^ 



(29) 



whose associated second order linear equation is a Riemann equation. More precisely, by 
Lemma 12.11 it can be written as a Legendre equation 

,2 



t^)u"{t) - 2tu\t) 



u{t) = 0, 



(30) 



with 



m 



2m 



and u being a solution of 



m 



fc2 



abk'^ 



Am? 



mc 



2abm'^ 



0. 



The difference of exponents in fl30|) is fi, /i and 2z/ — 1 and, therefore, we are under the 
hypotheses of Kimura's Theorem (see Appendix IB.1.1|) . 

Proposition 3.2. Legendre equation fl30|) is integrahle if and only if, either 

(1) ^^v^'Lorv^'L, or 

(2) ±fi, ±fi, ±(2i/ + 1) belong to one of the following seven families 



Case 


fi E 


u e 




(a) 


Z + i 


c 




ib) 




iZ ± i 


z + i 


(c) 


z±l 




Z + i 


id) 






^ + To 


(e) 


Z± i 

5 




^ + To 


(/) 


z± 1 

5 


iZ ± i 


Z + i 
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Proof. In Kimura's Theorem, the difference of exponents /i, /i and 2i/ + 1 correspond to the 
possibihties hsted above. Indeed, they are cases (i), (m.3), (ii.l2), (iil3) and 

(m.15) of Kimura's table (see Appendix IB.l.ip . For the case (i) we have 

• /i + yU + 2z/+lG2Z+l^// + z/GZ, 

• -fi + fi + 2u + le2Z+l^ueZ, 

• fi-fi + 2u + ie2Z + i^ueZ, 

• fi + fi- 2u-ie2Z + i^fi-ueZ. 

The rest of the cases can be proven in a similar way. 
(ii.l) We see that ±/i G | + Z and ±(2z/ + 1) G C and therefore /i G Z + | and u E C 
(ii.3) We consider that ±/i = / + |, ±/i = m + |, ±(2z/ + 1) = g + |, being l,m,q E Z. Take 

for instance /i G Z ± | and u E |Z ± |. Furthermore, / + m + g must be even and in 

consequence we obtain fi + u E Z + ^. 

(11.11) We have that ±/i = / + |, ±/i = m + |, ±(2z/ + 1) = g + | with l,m,q E Z. For 
example we consider /i G Z ± | and u E |Z ± |. Moreover / + m + g must be even 
and therefore we have that fi + u E Z + ^ . 

(11.12) Now we consider ±/i = / + |, ±/i = m + |, ±(2i/ + 1) = g + |, being l,m,q E Z. We 
take for instance /i G Z ± | and z/ G |Z ± |. Additionally / + m + g must be even and 
so fi + u E Z + j^. 

(11.13) Let be ±/i = / + |, ±/i = m + |, ±(2z/ + 1) = g + | with l,m,q E Z. Consider for 
example /i G Z± | and z/ G |Z± |. Furthermore, /-t-m + g must be even and therefore 
^ + u EZ + j^. 

(ii.l5) Consider ±/i = / + |, ±/i = m + |, ±(2z/ + 1) = g + |, being l,m,q E Z, take for 
instance /i G Z ± | and G |Z ± |. Moreover, / + m + g must be even and in 
consequence + u E Z + ^. 

Finally, observe that the difference exponents n, n and 2z/ + 1 do not satisfy the conditions 
(ii.2), (m.4), (m.5), (m.6), (m.7), (m.8), (m.9), (ii.lO) and (m.14) □ 



Now we deal with an example from [TTl [T2] . We consider the equation 

= A(x) + B(x)w, (31) 

dw 

where x is the dependent variable and w is the independent one. Now, for i? ^ 0, by means 
of the change of variable 

A 

and changing (ty, x) — > {x, y) (that is, we consider now x as the independent variable and y 
as the dependent one) equation (l3T|l is transformed into the Lienard equation 

dy I d f A\ , , 

^d^ = 5 + d^UJ'' ^''^ 

for any functions A and B ^ 0. In particular, for 

A = A{x) = a + bx + cx^, B = B{x) = a + f3x + ■yx'^, 
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equation (IHTj) falls into the Riccatti form 

dx 

—— = (a + aw) + (b + (3w)x + (c + 'yw)x'^. 
aw 

By Proposition 12. applying the transformation T it follows the reduced Riccatti equation 
(I13p and through the transformation TZ the normalized second order differential equation 
C," = p{x)$, with 

, , - 4a7 2 2a7 + 2ac - b(3 Aac - 6^ 67 - /3c 87^ 

p(x) = : X X 1 z r + — r-. (33) 

^ 4 2 4 2(7x + c) 4(7x + c)2 ^ ' 

Introducing the change r = 7X + c we get 

e = P(r)e, (34) 



where 



13^ - 4^7 2 207^ - 2ac7 - l3b-f + l3^c 6^7^ - 26/3c7 + l3^c^ 
^ V 4^2 

67 - f3c 372 



2r 4r2 ' 



and performing z = ■y ^—^^^t we arrive to 

r = ^{z)^, </.(z) = ^2 + 5i^ + |-52 + ^ + ^, (35) 

with 5i being algebraic functions in a, 6, c, a, /3 and 7. This equation (135 p is exactly the 
biconfluent Hewn equation whose integrability is analysed in Appendix B.2. 

Assuming /3 = 7 = we obtain a Lienard equation which is transformable into a reduced 
second order differential equation with r G C[a:] and deg(r) = 1. This means that the equation 
is not integrable (see [30] and section l3^ . As a particular case, we have a Lienard equation 
that can be reduced to the Riccati equation given in [?, equation 1.3.2.1], 

2yy' = {ax + b)y + 1. (36) 

Moreover, under some restrictions over the parameters of the biconfluent Heun equation 
one can obtain the Whittaker equation. For instance, the Lienard equation 

y^ = (A + biix) y+ {ex + b2ox'^) , 
falls into a Whittaker equation for some special values of the parameters. 

Remark 3.3. We would like to stress the following facts. 

(a) It is well-known that via the change z{x) = J f{x)dx (with inverse x = x{z)), the 
Lienard equation ( [F7| ) can be led into the equation 

yy'{z) = y + h{z), (37) 

with 

WW) 
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In a similar way, the change z{x) = J g{x)dx reduces (127|) to 

yy'{z) = h{z)y + 1, (38) 

with 

However, in general, these transformations do not preserve the differential field of the 
coefficients. This is why we do not use them to reduce equation (127|) . 
(b) Sometimes equation [27\ ) is called Abel equation of second kind since through the 
change y = l/w it is reduced to an Abel equation 

w' = —f{x)w^ — g{x)w^. 

In Cheb-Terrab [HI [12] the authors consider some families of Abel equations which are 
reducible to Riccati equations. 

3.4. Other families. Here we consider some special cases of Riccati equations. 

a) Polynomial Riccati equations 
The Riccati equation 

w' = r(x) — w'^, r{x) G C[x] 

has been studied by several authors (see, for instance [H EZJ [S3]). The Galois group of its 
associated second order linear differential equation is connected and can be either SL(2, C) 
or the Borel group [3l [2] . In the first case the tangent field associated to its Riccati equation 
has no invariant curves. In the second one, there is no rational first integral for its vector 
field. As an example, the reduced form for the triconfluent Heun equation is of this type and 
is given by 

e = P(x)e, p{x) = — + -62x' -6,x + ^- 5o. (39) 

b) Equations with finite Galois group 
Consider the polynomial Riccati vector field 

i = -x\x-lf 

y = a{x — ly + bx'^ + cx{x — 1) + x'^{x — l^v'^. ^ ' 

Its foliation was studied in |55] and is related to the following differential equation 

^ \^ 4x2 ^4(x-l)2^ 4x(x-l) ^ ' 

where A, /i and v are the differences of the exponents at 0, 1 and oo of the hypergeometric 
equation, see also Appendix [B] When the equation (HTj) is integrable their solutions are 
Legendre functions. For (A,/i, u) = (i, |, ^), the solutions of equation (HTj) are given by 



^1 = ^x2 -X (2x-l + 2Vx^ 
^2 = ^x^-x (2x-l + 2^x2" 



X 



X 
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and its differential Galois group is the Dihedral Group (case 2 of Kovacic algorithm). One 
can also obtain the tetrahedral, octahedral and icosahedral groups for (A,/i, z/) = (|,^,|), 
(A,/i, z/) = (|, |, |) and (A,/i, z/) = (|, |, |), respectively (case 3 of Kovacic algorithm). 
Related to fllT]) we have the equation (see [53] ) 



7x 



36z/2 



6x(x 



lAAx{x - 1) 



2/ = 0, 



(42) 



where the differential Galois group can be tetrahedral, octahedral and icosahedral depending 
on the values of u (indeed, for 1/3, 1/4 and 1/5 respectively). Even though explicit solutions 
of equations fHTl) and fH2l) are difficult to get in general, it can be proved the existence of 
rational first integrals for the associated vector field for values of the parameters z^ ^ {|, |, |}- 
For instance, solutions of equation for z/ = ^ (tetrahedral group) are 



„ _ 4/ 3/Z , -I sl 2V^+^+l-V3^^v/3 



In the case z^ = i (octahedral group) we obtain the solutions yi = ^^j"^ and y2 = /\ where 



72/3 



/i 

/2 
/3 



.6/ 3/ yi-x'+i 



ti+l+V3 



^+V3 



V3 



/T^-l 



/T=x + 1 



2/3 



+V3+2 



hi \ 2/3 



Finally, setting u = \ (icosahedral group), we get the solutions yi = '^^ and 



2/2 
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where cui = uJi{x) and a;2 = ^2(2;) are roots of the polynomial ^^afc(x)a;'^ and 



ai(x) 
a2(x) 
a^{x) 
a^i^x) 
a^{x) 



k=0 

102400- 11264x- llx^ 
-3686400X + 3679200x2 + 7200x^ 
479001600x2 - 476863200x^ - 2138400x^ 

30412800000x2 - 60445440000x^ + 29652480000x^ + 380160000x^ 
-821145600000x^ + 1597384800000x^ - 731332800000x^ - 44906400000x^ 
-3695155200000x^ + 11085465600000x'^ + 3695155200000x^ - 
11085465600000x^ 
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ae{x) = -492687360000000x^ + +1693612800000000x^ - 2124714240000000x'^ + 

1139339520000000x^ - 215550720000000x^ 
arix) = 8868372480000000x^ - 35473489920000000x^ + 53210234880000000a;^ - 

35473489920000000a;^ + 8868372480000000a;^ 
as{x) = -249422976000000000x*^ + 997691904000000000x^ - 14965378560000000000;^ + 

997691904000000000x^ - 249422976000000000x^° 
ag{x) = -4434186240000000000x'^ + 22170931200000000000x^ - 

44341862400000000000x^ + 44341862400000000000x^ - 

22170931200000000000x^° + 4434186240000000000x" 
aio(x) = 39907676160000000000x^ - 199538380800000000000X* + 

399076761600000000000a;^ - 399076761600000000000x^° + 

199538380800000000000a;" - 39907676160000000000x^^ 
au{x) = 

ai2(x) = 2176782336000000000000x^ - 13060694016000000000000^*^ + 

32651735040000000000000x^° - 43535646720000000000000a;" + 
32651735040000000000000x^2 - 13060694016000000000000a;^^ + 
2176782336000000000000x1^ 

c) Lame families. 

Let us consider in the Lame equation 

(Py_ ^ fix) dy _ n{n + l)x + B ^ ^ ^ ^^^^ 
dx"^ 2/(x) dx /(x) 

the case where /(x) = 4x^ — g2X — and parameters n, B, g2, gz such that the discriminant 
of /, namely 27(7| — f?!, is non-zero (see Appendix |B]) . Performing the transformation 

V = (44) 
cy 

with c = c(x) being any non-zero arbitrary rational function, it follows the family of Riccati 
equations associated to the Lame equation, 

, n{n + l)x + B ( f'{x) d , ^ 



For the Lame case, (i.l) of lB.3l with B = Bi, and from Remark lB.5l it turns out the existence 
of a polynomial solution of equation despite of the general solution of this equation is not 
algebraic. Furthermore, it is clear that for a fixed n we have associated Riccati equations f l45p 
with arbitrary degree, because the non-zero rational function c = c(x) at the transformation 
above is arbitrary. The Lame functions correspond here to algebraic solutions of (H5|) and 
moving n along the integers, we obtain algebraic invariants curves 

E'(x) 
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of unbounded degree. 

Hence, we have prove the following result. 

Theorem 3.4. For any fixed degree in the Riccati family [J^ ) associated to the Lame equation 
there exist invariant algebraic curves of any tangent field X of the corresponding foliation 
with unbounded degree. Furthermore, the first integral of X is not rational. 

Appendix 

Appendix A. Kovacic Algorithm 

This algorithm is devoted to solve the reduced linear differential equation (RLDE) ^" = pC, 
and is based on the algebraic subgroups of SL(2, C). For more details see [20]. Although 
improvements for this algorithm are given in 132] , we follow the original version given by 
Kovacic in |3U] . 

Theorem A.l. Let G be an algebraic subgroup o/SL(2, C). Then one of the following four 
cases can occur. 

(1) G is triangularizable. 

(2) G is conjugate to a subgroup of infinite dihedral group (also called meta-abelian group) 
and case 1 does not hold. 

(3) Up to conjugation G is one of the following finite groups: Tetrahedral group, Octahe- 
dral group or Icosahedral group, and cases 1 and 2 do not hold. 

(4) G = ^h{2X). 

Each case in Kovacic algorithm is related to each one of the algebraic subgroups of SL(2, C) 
and its associated Riccatti equation 

e' = r-e^ = {y/^-e){y/? + e), e = ^,r = p. 

According to Theorem lA.il we obtain four cases. Only for cases 1, 2 and 3 one can solve the 
differential equation RLDE and in case 4 one has not Liouvillian solutions for it. Kovacic 
algorithm can possibly provide one solution (^i), so the second one (^2) can be got through 

6 = / |. (46) 

Notations. For the RLDE given by 
we use: 

(1) Denote by T' be the set of (finite) poles of r, F' = {c G C : t(c) = 0}. 

(2) Denote by F = F' U {00}. 

(3) By the order of r at c G F', o(rc), we mean the multiplicity of c as a pole of r. 

(4) By the order of r at 00, o (roo) , we mean the order of 00 as a zero of r. That is 
o (roo) = deg(t) - deg{s). 
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A.l. The four cases. Case 1. In this case [a/t]^ and [a/^q^ means the Laurent series of A/r- 
at c and the Laurent series of a/t at oo respectively. Furthermore, we define e{p) as follows: 
if p G r, then e (p) G {+,—}• Finally, the complex numbers will be defined 

in the first step. If the differential equation has not poles it only can fall in this case. 

Step 1. Search for each c G F' and for oo the corresponding situation as follows: 

0. 
1. 

■ , then 



(Co) If o (r,) 


= 0, 


then 




0, 




(ci) If o(rJ 


= 1, 


then 




0, 




(C2) If o (rj 


= 2, 


and 


r = ■ ■ ■ + b{x - 


-c) 





in n ± i±vrT46 



c 



(cs) If o (rc) = 2v > 4, and 

r = (a (x - c)"'' + ... + d (x - + b{x - c)-(''+i) + ■ ■ ■ , then 

[a/t]^ = a (x - c)"" + ... + d{x - c)~^ , (yf = ^ + "W^ • 

(001) If o (roo) > 2, then 

(002) If o (roo) = 2, and r = ■ ■ ■ + + ■ ■ ■ , then 

[v^Joo = 0' «oo = ^ ■ 

(003) If o (roo) = — 2t' < 0, and 

r = {ax" + ... + df + te''"^ H , then 

[v^oo = + ••• + ^' " ^ (^^ ~ ■ 

Step 2. Find D 7^ defined by 



I m, G Z+ : m = - ), V (e (p))^,^ 1 

I cer' J 



If D = 0, then we should start with the case 2. Now, if j^D > 0, then for each m G -D we 
search u G C(x) such that 



eioo)[V^^ + J2 (c) [v^ , + «^^^)(x - c)-i) . 



' L • J 00 



Step 3. For each m G -D, search for a monic polynomial Pm of degree m with 

+ 2wP;, + {00' + - r)P„ = 0. 
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If one successes then = Pm^^^ is a solution of the differential equation RLDE. Else, 
Case 1 cannot hold. 

Case 2. Search for each c G F' and for oo the corresponding situation as follows. 

Step 1. Search for each c G F' and oo the sets ^ ^ and E^o 7^ 0- For each c G F' and 
for oo we define Ec G Z and E^o G Z as follows: 

(ci) If o (r,) = 1, then E, = {4}. 

(C2) If o (rc) = 2, and r = ■ ■ ■ + b{x — c)~^ + ■ ■ ■ , then 

= |2 + kVl + ^b ■.k = 0, ±2| . 

(cs) If o (rc) = V > 2, then Ec = {v}. 

(001) If o (roo) > 2, then = {0, 2, 4}. 

(002) If o (roo) = 2, and r = ■ ■ ■ + te^ + ■ ■ ■ , then 

= {2 + kVTT4b ■.k = 0, ±2| . 

(003) If o (roo) = V <2, then E^o = {v}. 
Step 2. Find D ^ defined by 

Z}=<|mGZ+: m = ]- icoo - | , Vcp G -Ep, p G F 



cer' 



If D = 0, then we should start the case 3. Now, if j^D > 0, then for each m G -D we search 
a rational function 6 defined by 



9 



-E- 



Step 3. For each m ^ D, search a monic polynomial Pm of degree m, such that 

p^' + s^p;;; + (3^' + 39"^ - 4r)p^ + {9" + m' + 9^ - Ar9 - 2/) p^ = 0. 

If Pm does not exist, then Case 2 cannot hold. If such a polynomial is found, set = 9 + P' / P 



and let a; be a solution of 



+ 0a; + ^ (0' + 02 - 2r) = 0. 



Then 6 = e^"" is a solution of the differential equation RLDE. 

Case 3. Search for each c G F' and for oo the corresponding situation as follows: 

Step 1. Search for each c G F' and oo the sets -Ec 7^ and E^ ^ 0. For each c G F' and 
for 00 we define E^. GT, and E^o C Z as follows: 

(ci) Ifo(rc) = l, thenPc = {12}. 

(C2) If o (rc) = 2, and r = ■ ■ ■ + h{x — c)"^ + ■ ■ ■ , then 

Ec = je + k^l + Ah : k = 0, ±1, ±2, ±3, ±4, ±5, ±6| . 
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(oo) If o (roo) = f > 2, and r = ■ ■ ■ + bx"^ + ■ ■ ■ , then 

f 12k ^ 

= j6 + — Vl + 46: A; = 0,±1,±2,±3,±4,±5,±6 ^ , nG {4,6,12}. 



Step 2. Find D 7^ defined by 



D = { m ^ Z+ : m = — I e, 



00 



^Cc^ ,Vep e -Ep, p G r| . 



In this case we start with = 4 to obtain the solution, afterwards n = Q and finally n = 12. 
If -D = 0, then the differential equation has not Liouvillian solution because it falls in the 
case 4. Now, if > 0, then for each m & D with its respective n, search a rational function 

cer' 

and a polynomial 5* defined as 

S=l[{x-c). 

cer' 

Step 3. Search for each m & D, with its respective n, a monic polynomial Pm = P of 
degree m, such that its coefficients can be determined recursively by 

P_l = 0, Pn = -P, 



P,_i = -SP: - {{n -z)S'- S9) Pi-{n- i) {i + 1) ^VP,+i, 

where i G {0, l...,n — l,n}. If P does not exist, then the differential equation has not 
Liouvillian solution because it falls in Case 4. Now, if P exists search u such that 

^(n-z)! 

then a solution of the differential equation the RLDE is given by 
where u is solution of the previous polynomial of degree n. 



Appendix B. Some Special Functions 
B.l. Hypergeometric families. 
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B.1.1. Kimura's Theorem. The hypergeometric (or Riemann) equation is the more general 
second order hnear differential equation over the Riemann sphere with three regular singular 
singularities. If we place the singularities at x = 0, 1, oo it is given by 

cPy f \ — a — a' 1 — j — dy 



dx"^ \ X X — 1 J dx 

' aa' 77' /3/3' — aa'^^' 



(47) 



x"^ (x — 1)^ x{x — 1) 

where (a, a'), (7,7'), {P^P') are the exponents at the singular points and must satisfy the 
Fuchs relation a + a' + 7 + 7' + /3 + /5' = 1. 

Now, we will briefly describe Kimura's Theorem that provides necessary and sufficient 
conditions for the integrability of the hypergeometric equation. Let be A = a — a', /i = /3 — 
and z/ = 7 — 7'. 

Theorem B.l (Kimura, [SH])- The hypergeometric equation ( (7^ is integrable if and only if 
either 

(i) At least one of the four numbers A + yU + z/, —A + ^ + v, X — fi + u, X + fi — u is an 
odd integer, or 

(a) The numbers X or —A, /i or —fi and v or —v belong (in an arbitrary order) to some 
of the following fifteen families 



1 


1/2 + / 


1/2 + m 


arbitrary complex number 




2 


1/2 + / 


1/3 + m 


1/3 + q 




3 


2/3 + / 


1/3 + m 


1/3 + q 


I + m + q even 


4 


1/2 + / 


1/3 + m 


1/4: + q 




5 


2/3 + / 


1/4 + m 


1/4: + q 


I + m + q even 


6 


1/2 + / 


1/3 + m 


1/5 + g 




7 


2/5 + / 


1/3 + m 


1/3 + g 


I + m + q even 


8 


2/3 + / 


1/5 + m 


1/5 + g 


I + m + q even 


9 


1/2 + / 


2/5 + m 


1/5 + g 


I + m + q even 


10 


3/5 + / 


1/3 + m 


1/5 + g 


I + m + q even 


11 


2/5 + / 


2/5 + m 


2/5 + g 


I + m + q even 


12 


2/3 + / 


1/3 + m 


1/5 + g 


I + m + q even 


13 


4/5 + / 


1/5 + m 


1/5 + g 


I + m + q even 


14 


1/2 + / 


2/5 + m 


1/3 + q 


I + m + q even 


15 


3/5 + / 


2/5 + m 


1/3 + q 


I + m + q even 



Here l,m and q are integers. 

B.l. 2. Confluent hypergeometric. The confluent Hypergeometric equation is a degenerate form 
of the Hypergeometric differential equation where two of the three regular singularities merge 
into an irregular singularity. The following are two classical forms. 

• Rummer's form 

y"+^^y'--y = 0, a, C E C (48) 
X X 
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Whittaker's form 



y 



1 K V - 1 
4 X 



(49) 



where the parameters of the two equations are hnked by and yU = f Furthermore, 



using the expression ffTOj) . we can see that the Whittaker's equation is the reduced form of the 
Kummer's equation fHSjl . The Galoisian structure of these equations has been deeply studied 
in I1I1I22]. 

Theorem B.2 (Martinet & Ramis, [H]). The Whittaker's differential equation (l49l) is in- 
tegrable if and only if either, k + n ^ \ + or k — fi E | + or —k + /i G | + N, or 

-K - fi e I + N. 

The Bessel's equation is a particular case of the confluent Hypergeometric equation and is 
given by 



y +-y 

X 



-y = 0. 



(50) 



Under a suitable transformation, the reduced form of the Bessel's equation is a particular 
case of the Whittaker's equation fH9|) . 



Corollary B.3. The Bessel's differential equation fl50p is integrable if and only if n G | + 

B.2. Heun's families. The Heun's equation is the generic differential equation with four 
regular singular points at 0, 1, c and 00. In its reduced form, the Heun's equation is y" = 
r{x)y, where 



r{x) 



A 

X 



B 



+ 



C 



^ a/3 a'j 5rjh 



X — 1 X — c 
B ■ 



D 

~2 



E 



+ 



X 



(X-1)2 (X-C)2' 

/37 57]{h-l) 



(51) 



2{c-r 



c- 1 



(3'y 6ri{c — h) ^ a /a 
2c ' 2(c- 1) ~ c(c- 1) ' ~ 2 V 2 



2 V 2 



F = |(|-l), with a + /3 + -f-6-r] = l. 
To our purposes we write the determinant 11^4.1(0, b, u, v, C,, w) as in [22] : 



w u 

dS,w + l V 2{u + b) 

(rf-l)C w + 2{v + a) 3(u + 25) 

(rf-2)C w + 3(i; + 2a) 4(u + 36) 



2C w+(d-l)(w + (d-2)a) d(u+(d-l)6) 
^ w + + (rf- l)a) 
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B.2.1. Biconfluent Heun. The equation 

(52) 

is the well known biconfluent Heun equation which has been deeply analyzed by Duval and 
Loday-Richaud in [22l p. 236]. 

Theorem B.4. [22j. The biconfluent Heun equation fl52l) has Liouvillian solutions if and 
only if it falls in Case 1 of Kovacic algorithm and one of the following conditions is fulfilled. 

(1) 5l = l, 5:i = Q and ^2 G 2Z + 1. 

(2) 5o = 1; ^3 7^ and 82 G 2Z* + 1 with \52\ > 3, and if e = sign52, then 

n(|52|-i)/2 1, 2, e6i, -2e, e6i - =0. 

(3) 60 ^ ±1, ±60 ±62 e 2Z* and if Eq^e^ G {±1} are such that E00S2 - SqSo = 2(i* G 2W 
then 

n^. ^0, 1,1 + £0^0, ^oo^i, -2£oo, ^(£00^1(1 + £oSo) - 63)^ = 0. 

B.3. Lame equation. The algebraic form of the Lame Equation is [171 156] 

d^y f'{x) dy n{n + l)x + B ^ 

dx^^2f{x)dx fix) ^ ' ^ ^ 

where f{x) = — g2X — g^, with n, B, g2 and g^ parameters such that the discriminant 
of /, A = 27 gl — gl, is non-zero and, therefore, it has no multiple roots. This equation is a 
Fuchsian differential equation with four singular points over the Riemann sphere: the roots 
ei, 62, 63 of / and the point of the infinity. 

The mutually-exclusive known cases of solutions in closed form of the Lame equation (153|) 
are the following. 

(i) The Lame-Hermite case (see [771156] ). We have n G N and arbitrary parameters -B, g2 
and g^,. 

(a) The Brioschi-Halphen-Crawford case (see [23|H7]). We have n G N such that m := n+ 
i G N and parameters B, g2 and g^ satisfying an algebraic condition Qm{.g2/^, 9zl^-, B) = 
0, where Qm G Z[g2/4:, g^/A, B] is a poljTiomial of degree m in 5, known as the Brioschi 
determinant. 

(Hi) The Baldassarri case (see [1]). One asks n to satisfy that n + | G |Z U U |Z — Z 
besides some additional (involved) algebraic restrictions on the other parameters. 

It is possible to prove that the only integrable cases of the Lame equation are cases {i)-{iii) 
above, see [12]. In (ii) and (Hi) the general solution of (153|) is algebraic and the Galois group 
is finite. Case (i) splits in the following two subcases [13 [56]. 

(z.l) The Lame case. For a fixed integer n, this equation admits a solution (called Lame 
function) of the form 

3 

E{x) = l[{x-e,)'^Pm{x), (54) 

i=l 
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being Pm. a monic polynomial of degree m = n/2 — {ki + k2 + ks) and ki G {O, |}, 
i = 1,2, 3. Since m G N, eight different possibilities regarding n appear: If n is even, 
we have ki = k2 = k^ = oi just one zero kf, if n is odd, we could have all non-zero 
/cj's or combinations with exactly one non-zero fcj. All these possibilities give rise to 
classes of Lame functions. Concerning parameter B, it must be one of the m + 1 
different roots Bi, . . . , B^+i of certain irreducible polynomial of degree m + 1, with 
all its roots real and simple [23]. Furthermore, the numbers Bi are reals. 
(i.2) The Hermite case. Here we are not in case (i.l) and n is an arbitrary natural number. 
We are also fall in case 1 of Kovacic algorithm, but with a diagonal Galois group. 

Remark B.5. We notice that the polynomial Pm in (i.l) satisfies a second order linear 
differential equation similar to the one that appears in the first case of Kovacic algorithm. In 
fact it is possible to obtain the above passing into normal form and applying Kovacic algorithm. 
Therefore the second linear independent solution is not algebraic and the associated Riccati 
equation has no rational first integral. 
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